可压缩多介质流动问题的高精度数值模拟方法

在可压缩流动问题的数值模拟领域,激波的高分辨率计算已取得重要进展。但是包含物质界面的可压缩多介质流动的数值模拟还存在诸多数值挑战,主要表现为界面处数值耗散过大和非物理振荡等问题。界面处流体性质的不连续性是造成可压缩多介质流动问题物理建模与数值方法困难的主要原因。为了建立一套高效的可压缩多介质流动问题的高精度数值模拟方案,本文从数值框架的选择、非守恒方程相容离散、高精度有界格式构造、界面压缩、界面-单介质分区计算方法等多个维度出发,综述近几年我们在可压缩多介质流动问题高精度数值模拟方法方面的研究进展。通过上述多个维度的工作,我们建立了一套适用于可压缩多介质流动问题的低耗散、基本无数值振荡的高精度欧拉数值方法,并成功应用于可压缩多介质大变形流动和界面不稳定性诱导湍流混合等问题的数值模拟。相关数值方法研究成果已集成至武器物理等领域工程数值模拟软件中,为相关工程任务提供了重要技术支撑。     

半球谐振子疲劳寿命分析与仿真

半球谐振子的寿命和损伤是直接影响高精度半球谐振陀螺使用时长和安全性的重要因素。目前国内加工的半球谐振子所用的熔融石英玻璃材料主要依靠进口，采用传统的疲劳寿命实验确定方法成本过于昂贵，因此需要利用疲劳分析软件对半球谐振子的疲劳寿命进行分析。文章通过半球谐振陀螺应力分析，采用ANSYS软件对熔融石英半球谐振子进行应力分析仿真，确定因残余应力所引起的疲劳危险部位，并在疲劳部位进行裂纹扩展分析，得到 、 和 型应力强度因子，以使在半球谐振子结构设计和使用过程中对易疲劳的部位进行有效监控和预防。     

基于DOA聚类算法的ARM抗有源诱偏技术研究

有源诱偏技术是雷达对抗反辐射导弹的一种重要手段,可以大大降低反辐射导弹的作战效能。为提高对敌目标的精确打击能力,在分析闪烁诱饵诱偏原理的基础上,针对有源诱偏干扰下被动雷达测角精度与稳定度不高的问题,通过对有源诱偏信号的时域特征进行分析,提出一种基于脉冲前沿检测的DOA聚类分选算法,找出前沿超前的辐射源信号,实现了高性能抗有源诱偏干扰和高精度抗干扰测向。仿真结果表明,该技术可以有效对抗四点源有源诱偏干扰。     

航天嵌入式软件数组越界缺陷特征研究

根据统计,数组越界是航天嵌入式软件开发过程中出现最多且最容易被遗漏的缺陷类型之一.目前自动化检测数组越界多基于抽象解释、符号执行、程序模型检验等方法,这些方法在误报、漏报、可扩展性等方面的表现依赖于软件及缺陷特征.分析了近三年航天嵌入式软件第三方测试中发现的94个数组越界问题,从缺陷模式和缺陷表现形式两方面分析得出10项航天嵌入式软件数组越界缺陷特征,并提出对设计具体检测方法关键的若干启示.进一步基于这些特征和启示探讨了数组越界检测算法针对中断驱动型程序的改进方向.     

一种基于OPTICS聚类的仅测角导航目标检测算法

仅测角自主导航方法具有设备简单、复杂度低，功耗低的优点，在空间任务中具有广泛的应用前景.针对中远距离下空间目标特征少的特点，提出了一种利用基于OPTICS聚类算法的空间目标检测方法，可用于仅测角导航过程中的目标检测.对原始星图进行预处理提取星点及目标点，并结合星图识别的结果选择部分帧，使用经过改进的OPTICS聚类方法获得目标运动轨迹.最后，使用本文中的算法对软件仿真出的含有目标的高精度星图进行处理验证了算法的可行性.在卫星相对于空间目标抵近过程中，目标检测的水平误差及垂直误差小于0.15°的帧数分别占到了85.4%以及99.6%.相比AVANTI实验中的目标检测方法，减少了在轨任务中相关参数的调节，进一步提升了算法的自主性.     

微纳聚合体卫星的对称式重构规划算法

为解决微纳聚合体卫星变构过程中，聚合体各部分之间的动力学耦合导致的卫星整体姿态翻转或紊乱，提出了一种对称式变构规划算法。首先建立了铰链约束下的漂浮基多刚体系统动力学模型，研究了重构过程中各运动模块对本体姿态的影响，提出当每两个运动模块位置与转动方向满足对称性条件时，两者对本体姿态的影响可在一定程度上相互抵消。基于上述理论，在A*算法中引入最优分配度量和对称性判定，设计了并行对称重构规划算法。仿真结果表明，该规划算法可实现重构过程中多模块并行、对称运动，重构过程总步数较少，运动模块对本体姿态的影响小。     

短螺旋型燃烧室旋流结构及气动边界发展规律数值模拟研究

短螺旋型燃烧室的头部轴线沿发动机周向与发动机轴线偏转一定角度,可以在保证燃烧效率不变的条件下有效缩短燃烧室轴向长度。短螺旋型燃烧室流场的最大特点是旋流流动单侧受限。为了研究头部安装角α变化对燃烧室旋流流动特性的影响,基于数值方法对短螺旋型燃烧室进行了计算分析。结果表明:随着α变化,旋流器下游旋涡依次出现对称环状、马蹄状、环状结构;随着α增加,气动边界逐渐出现并抑制旋流的周向扩张,导致流场出现不同的旋涡结构;不同α下的切向角动量随轴向距离增加而衰减,但α为35°和45°时,气动边界在非受限侧出现并对旋流产生约束,角动量衰减变慢;当α为0°,15°,35°,45°时,燃油液滴依次集中分布在旋流器下游两旋涡边缘、侧壁面和头部端壁、非受限旋涡边缘。本文研究了不同α下的旋涡结构及气动边界沿轴向的演化过程,为短螺旋型燃烧室进一步的设计与优化提供基础。     

A comprehensive survey on the methods of angle of attack measurement and estimation in UAVs

Angle of Attack (AOA) is a crucial parameter which directly affects the aerodynamic forces of an aircraft. The measurement of AOA is required to ensure a safe flight within its designed flight envelop. This paper intends to summarise a comprehensive survey on the measurement techniques and estimation methods for AOA, specifically in Unmanned Aerial Vehicle (UAV) applications. In the case of UAVs, weight constraint plays a major role as far as sensor suites are concerned. This results in selecting a suitable estimation method to extract AOA using the available data from the autopilot. The most feasible and widely employed AOA measurement technique is by using the Multi-Hole Probes (MHPs). The MHP measures the AOA regarding the pressure variations between the ports. Due to the importance of MHP in AOA measurement, the calibration methods for the MHP are also included in this paper. This paper discusses the AOA measurement using virtual AOA sensors, their importance and the operation.     

Unsteady wakes-secondary flow interactions in a high-lift low-pressure turbine cascade

Detailed experimental measurements were conducted to study the interactions between incoming wakes and endwall secondary flow in a high-lift Low-Pressure Turbine (LPT) cascade. All of the measurements were conducted in both the presence and absence of incoming wakes, and numerical analysis was performed to elucidate the flow mechanism. With increasing Reynolds number, the influence of the incoming wakes on suppressing the secondary flow gradually increased owing to the greater influence of incoming wakes on reducing the negative incidence angle at higher Reynolds numbers, leading to a lower blade loading near the leading edge and suppression of the Pressure Side (PS) leg of the horseshoe vortex. However, the effect of unsteady wakes on suppressing the profile losses gradually became weaker owing to the reduced size of the Suction Side (SS) separation bubble and increased mixing loss in the free-flow region at high Reynolds numbers. Incoming wakes clearly improved the aerodynamic performance of the low-pressure turbine cascade at low Reynolds numbers of 25,000 and 50,000. In contrast, at the high Reynolds number of 100,000, the profile loss at the midspan and mass-averaged total losses downstream of the cascade were higher in the presence of wakes than in the absence of wakes, and the unsteady wakes exerted a negative influence on the aerodynamic performance of the LPT cascade.     

Effect of bristle pack position on the rotordynamic characteristics of brush-labyrinth seals at various operating conditions

Over the last few decades, the research on the effect of bristle pack position on the rotordynamic characteristics of the brush-labyrinth seals is not sufficient. To this end, two kinds of brush-labyrinth seals for the bristle pack element installed upstream of the labyrinth teeth named BSU and installed downstream of the labyrinth teeth called BSD were used to investigate the effect of bristle pack position on the rotordynamic characteristics of the brush-labyrinth seals. Using the numerical model combining the porous medium model and the whirling rotor method, the rotordynamic characteristics of the BSU and BSD at various operating conditions including four kinds of pressure ratios, five kinds of inlet preswirl speeds and four kinds of rotor spinning speeds were conducted. The obtained results show that the effects of operating conditions on rotordynamic coefficients for the different seal configurations are different. The direct stiffness, cross-coupled stiffness and direct damping of the BSU are lower than those of the BSD. The rotordynamic coefficients of the BSU are more insensitive to the operating conditions variation. From the perspective of the seal stability, the BSU is a better brush-labyrinth seal configuration at high pressure ratio, high positive preswirl or high rotor spinning speed conditions. While in the case of low pressure ratio, low positive preswirl or low rotor spinning speed conditions, the BSD is a better choice.